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We study the physical limitations placed on CMB temperature and polarization measurements of 
the initial power spectrum by geometric projection, acoustic physics, gravitational lensing and the 
joint fitting of cosmological parameters. Detailed information on the spectrum is greatly assisted 
by polarization information and localized to the acoustic regime k = 0.02 — 0.2 Mpc~^ with a 
fundamental resolution of Ak/k > 0.05. From this study we construct principal component based 
statistics, which are orthogonal to cosmological parameters including the initial amplitude and 
tilt of the spectrum, that best probe deviations from scale-free initial conditions. These statistics 
resemble Fourier modes confined to the acoustic regime and ultimately can yield ~ 50 independent 
measurements of the power spectrum features to percent level precision. They are straightforwardly 
related to more traditional parameterizations such as the the running of the tilt and in the future 
can provide many statistically independent measurements of it for consistency checks. Though we 
mainly consider physical limitations on the measurements, these techniques can be readily adapted 
to include instrumental limitations or other sources of power spectrum information. 



I. INTRODUCTION 

High precision cosmic microwave background (CMB) 
temperature anisotropy measurements from WMAP have 
ushered in a new era in the determination of the initial 
conditions for structure formation and their implications 
for the physics of the early universe (ij . Hints of devia- 
tions from power law initial conditions from WMAP in 
the form of marginally significant coherent glitches over 
several multipoles and running of the tilt can be tested 
by the Planck satellite ^ and definitively measured with 
a high precision CMB polarization experiment such as 
the envisag ed Inflation Probe CMBPol % . 

Significant deviations would imply that the initial 
power spectrum is not scale free and hence traditional 
descriptions such as a constant running of the tilt, mo- 
tivated by the small deviations expected from slow-roll 
inflation (e.g. Q), can lead to misinterpretation and bi- 
ases in parameter determination. 

A large body of work already exists on more general 
approaches. They are typically based on flat bandpow- 
ers, discrete sampling with interpolation for continuity or 
smoothness, wavelet expansions or more generally some 
predetermined set of windows on, or effective regulariza- 
tion of, the initial power spectrum |E IE S 13 ■ 

In this paper, we study the fundamental physical lim- 
itations placed on the measurement of the initial power 
spectrum by geometric projection, acoustic physics, grav- 
itational lensing, and the joint fitting of cosmological pa- 
rameters. We then use this study to develop complemen- 
tary techniques that are optimized with respect to the 
physics of the CMB rather than prior assumptions about 
the form of the initial power spectrum. 

The outline of this paper is as follows. In we study 
the processes that transfer fluctuations from the initial 
curvature perturbation to the observed temperature and 



polarization flelds and quantify the limitations imposed 
by acoustic physics and geometric projection. In Wl\ we 
introduce a discrete parameterization of the initial spec- 
trum that is matched to the further limitations imposed 
by gravitational lensing. In i)IVI we employ Fisher ma- 
trix techniques to study the effect of joint estimation of 
the initial spectrum and cosmological parameters. In 
we construct modes which best constrain from scale-free 
initial conditions out of its principal components. We 
conclude in ijVIl 



II. 



TRANSFER FUNCTIONS AND 
PROJECTION 



Acoustic and gravitational physics generate CMB tem- 
perature fluctuations Q = AT{h)/T and polarization 
fluctuations E{ii) ^ in the angular direction fi from an 
initial spatial curvature fluctuation C(x), which we will 
specify in the comoving gauge. The observable angular 
power spectra are deflned by the two point function of 
their multipole moments 
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where A, A' £ 6, £^ and consequently are related to the 
initial curvature power spectrum 

(C(k)C(k')) = (2vr)35(k - k')^A2(fc) . (2) 

For slow-roll inflation the initial spectrum is given by 
(e.g. see j^) 
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where Hi is the Hubble parameter, e = 3{1 + Wi)/2, and 
Wi is the equation of state, all evaluated when the fc-mode 
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FIG. 1: Transfer functions in temperature and polarization 
Contours are spaced by factors of 4 with blue (red) denoting 
fluctuations with the same (opposite) sign as Note the pro- 
jection ridge along the diagonal, the out-of-phase temperature 
and polarization acoustic oscillations along the ridge begin- 
ning with the Sachs-Wolfe (SW) effect for the temperature, 
the reionization bump for the polarization, and the potential 
decay terms from the dark energy (ISW) and radiation for 
the temperature. The polarization transfer is cleaner than 
and complementary to the temperature transfer. 



exited the horizon during inflation. To the extent that 
Hi and Wi are constant, the spectrum is scale invariant; 
to the extent that the evolution is small across observable 
scales, the spectrum is a scale-free power law. 

The mappings between the two sets of power spectra 
are specified by transfer functions 



27r 



{k)Al{k) 



(4) 



Note that unlike the familiar transfer function of the mat- 
ter power spectrum, the CMB transfer functions are in- 
herently two-dimensional since t hey convert spatial fluc- 
tuations to angular fluctuations [W^ . 

These transfer functions are obtained by numeri- 
cally solving the linearized Einstein-Boltzmann equa- 
tions. The publically available CMBFast [13 and CAME 
|l2j | integral codes do return these functions in princi- 
ple. However they implicitly assume that the initial con- 
ditions are smooth so that the observable power spec- 
tra can be rapidly calculated by interpolation in £ and 
k. This interpolation can cause problems if the initial 
power spectrum contains features that are comparable 
to the gridding, e.g. typically A£ = 50. Traditional 
Boltzmann hierarchy codes 0,0,^3 on the other hand 
solve the equations for every £ < £max but typically sam- 
ple sparsely in the computational ly c ostly fc-modes and 
so require explicit smoothing in £ |l6l | . 

We have instead developed an independent, fast ver- 
sion of a Boltzmann hierarchy code whose speed with- 
out smoothing is comparable to the public codes with 
an evaluation at every £. Specifically the code solves 
the linearized equations |l7l | and the multilevel-atom cal- 
ibrated 2-level hydrogen and helium recombination equa- 
tions in the comoving gauge out to £niax = 6000 
with absorptive boundary conditions. We have checked 
that the agreement with CMBFast v4.3 with RECFast 
for smooth initial conditions is typically at the few 10^'^ 
level and comparable to the numerical error in CMBFast 

For this work, we oversample the transfer functions 
with 3750 logarithmically spaced fc-modes from k — 
20-4.35 _ 10-1-35 Mpc-i {Sk/k = 1.8 X 10-3) and 2500 fc- 
modes from 10-1-35 _xo-o.35 Mpc-^ {5k/k = 0.9 xlO^^). 
Figure H shows the transfer functions for the fiducial 
flat cosmology with dark energy density Hue — 0.72, 
dark energy equation of state w — —1, matter density 
^mh^ = 0.145, baryon density Vthh"^ — 0.024 and reion- 
ization optical depth r = 0.17. 

The most notable feature of the transfer functions is 
the main projection ridge along the diagonal at fc_D, — 
where the angular diameter distance to recombination 
-D, — 13.8Gpc in the fiducial model. It is instructive 
to consider first the simple Sachs- Wolfe (SW) limit 
of purely gravitational effects in a fiat universe. Here 
the temperature fiuctuation is a projection of 1/3 of the 
Newtonian gravitational potential or equivalently 1/5 of 
the comoving curvature. 
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FIG. 2: CMB angular and spatial power spectra for the fidu- 
cial cosmology and spectrum = (5.07x10^^)^. The spatial 
power spectra A^(fc) are plotted with k — l/Dt.. The low k 
cutoff is due to Imin ~ 2 and corresponds to superhorizon 
scale modes. Asymmetry in the projection leads to ^ < kD^ 
power that is much more pronounced in the temperature than 
the polarization. 



FIG. 3: Angular spectrum of an individual k mode for a 
trough fc,„~5/2 = 0.05 Mpc^^ and a peak fcm~4 = 0.082 
Mpc~^. At acoustic troughs, the temperature perturbations 
are dominated by the Doppler effect which is smeared out in 
angle by the projection; here the polarization is at its maxi- 
mum and yields a sharp projection of power. 



Terefore comparing Eqn. ^ Q and Q), we obtain 

5 

Tf{k) = 0. (6) 

By further assuming a fiducial model with a scale- 
invariant initial power spectrum 

^c(fc)|fid='5c'-(5.07xlO-Y, (7) 
we can evaluate the integral over In k to obtain 
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The sums over In £ yield the spatial power spectra 
Ae(fc) - E ]TP{k)Tpik)Al « ^ , » 1) , 

Alik) ^ E V(fc)Tf(fc)A2=0, (9) 

which quantifies the contribution per log interval to the 
field variance {X{h)X{h)) 0|. Thus geometric projec- 
tion represents a power conserving mapping described by 



the spherical Bessel function ji{kD^). Under the SW ap- 
proximation the transfer function of Fig. ^ would appear 
as a diagonal band that narrowed with increasing k. 

The actual transfer functions differ from the SW ap- 
proximation in interesting ways. Figure [21 compares the 
£-space and fc-space power spectra for the fiducial model. 
The oscillations at high k oi £ are the familiar acous- 
tic oscillations. For the temperature transfer function, 
the oscillations mainly trace the local plasma temper- 
ature at recombination and modulate the transfer as 
— cos(fcs*). Here is the sound horizon at recombi- 
nation; = 144.5 Mpc"-'^ in the fiducial model. We will 
number the mth acoustic temperature peak as having 
wavenumber km — rmr/s^. 

Figure 121 shows a mode near the km=i- Though the 
projection has a peak at £ k, kD^, with a narrow width 
of A£ 10 — 20, there is a strong asymmetry in the tails 
that skews the transfer of broadband power to £ < kD^, . 
This asymmetry reflects the fact that a wave that is 
nearly parallel to the line-of-sight will have a projected 
wavenumber £ <C kD,, as the ji(kD<,) factor in the trans- 
fer implies. The asymmetry in the projection places the 
angular peaks £m at values that are a few percent less 
than k,nD^ The asymmetry has a larger affect 

on the broadband power on small scales. Since photon 
diffusion during recombination cuts off the spectrum at 
k — ku w 0.1 Mpc"""^, a given I in the damping tail shows 
a deficit of broadband power from the oscillating tails of 
high k modes. 
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Even in fc-space, there is finite power at half integral 
m. Here the plasma velocity induces a Doppler effect that 
modulates the transfer function as — sin(A:s*). Moreover 
the ^-structure of the transfer function differs qualita- 
tively from the SW result. Since the velocity transverse 
to the line of sight yields no Doppler effect, the transfer 
function distributes fc-power as j'^ instead of jg |lOl | . In 
Fig. Owe show a mode with m = 5/2. Notice that there 
is no concentration of power near i = kD^. Since power 
in these modes is distributed only broadband in £, re- 
covery of features in the initial power spectrum at these 
wavenumbers from the temperature power spectrum will 
be severely limited. 

At fc < tt/s*, effects between recombination and the 
present further broaden the transfer of power. As the 
photons traverse decaying potentials wells (hills) during 
the end of radiation domination and the onset of dark en- 
ergy domination, they suffer a net gravitational blueshift 
(redshift). We call latter the integrated Sachs Wolfe ef- 
fect and the former the radiation effect (also known 
as the early ISW effect 0). Since the distance to the 
observer is now D < D^, the projection takes the power 
to lower £ or larger angles. Conversely a single i acquires 
power from a broad range in k. We show individual £- 
elements of the transfer function in Fig. 0] for £ = 2, 
£m=2 = 531 and £m=3 = 812. Note the broad distribution 
in In k for contributions to the temperature quadrupole 
£ = 2. Even at the mth peak in £ there are secondary 
contributions from the wavenumbers near fc„+2/2 due to 
the broad projection of the Doppler effect. 

The acoustic and projection limitations of the tem- 
perature power spectrum are largely removed with po- 
larization information. Near the end of recombination, 
photon diffusion creates local quadrupolar temperature 
anisotropy. Thomson scattering of radiation with a 
quadrupolar temperature anisotropy linearly polarizes it. 
The acoustic polarization transfer is therefore modulated 
according to velocity gradients or (fc/fc/^) sin(fcs*) and 
like the Doppler effect is out of phase with the local tem- 
perature. Unlike the Doppler effect, the projection of 
polarization is sharp (see Fig.|2Jl. The spherical decom- 
position of plane wave E'-polarization projects as 
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leading to a nearly diagonal transfer function (see Fig.^. 
Furthermore there is only one source of acoustic polariza- 
tion unlike the temperature and so the power spectrum 
in k space A|;(fc) possesses zeros at km, where m is an 
integer. These sharp features in k are only moderately 
smoothed by the projection. Finally the low k, low £ 
end of the polarization transfer function is dominated by 
rescattering during reionization. The projection remains 
fairly sharp in comparison to the ISW dominated tem- 
perature projection though the shorter distance to the 
reionization epoch leads to a slightly larger £ for a given 
k. 
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FIG. 4: Spatial spectrum of an individual £ mode for a scale 
invariant initial spectrum. The temperature quadrupole £ = 2 
receives ISW contributions across a wide range in k while the 
polarization quadrupole has a much tighter distribution. At 
the acoustic peaks {£m=2 ~ 531, £m=3 = 812) the temperature 
distribution is tighter but still receives contributions around 
~ km+1/2 due to the broad Doppler distribution. This effect is 
more prominent in the polarization where the peaks represent 
deep minima that are filled in by the projection broadened 
peaks. 



In summary, the temperature and polarization fields 
are complementary in the transfer function, both in the 
out-of-phase acoustic power in k and in the sharpness 
of the geometric projection to £. With both spectra and 
the cross correlation, there is high sensitivity to the initial 
power spectrum in the whole acoustic range k ~ 0.02—0.2 
Mpc^^ with projection limiting the sharpness of features 
in £ to A£ > 10 — 20. The reionization bump in the polar- 
ization also offers a cleaner probe of the long wavelength 
modes. Although it is still severely limited by cosmic 
variance and the determination of the reionization his- 
tory from smaller scale modes, polarization can in prin- 
ciple provide an incisive test for long wavelength power 
in our horizon volume. It can then determine whether 
the anomalously small low order temperature multipoles 
found by COBE and WMAP are due to a lack of long 
wavelength power in our horizon volume or a chance can- 
cellation with more local effects. 



III. DISCRETE REPRESENTATION AND 
GRAVITATIONAL LENSING 

Due to the geometric projection and, as we shall see, 
gravitational lensing, sharp features in the initial power 
spectrum will not be preserved in the observable angular 
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power spectra. For these reasons the transfer function 
cannot be inverted directly to obtain the initial power 
spectrum from the measured angular power spectra even 
with perfect measurements. Furthermore, as we shall see 
in the next section, the transfer function itself depends on 
cosmological parameters that must be jointly estimated 
from the data. 

Given the physical limitations of the transfer function 
it is convenient to represent the initial power spectrum by 
a discrete set of parameters. This discretization will also 
facilitate the treatment of gravitational lensing below. 
Specifically, we represent the initial power spectrum as 

\nAl{k)=ln^l + Y,P^W^{lIlk), (11) 

i 

where W is the triangle window 

/ I In fc — In fc, I \ 
VF,(lnfc) = maxh-|^^-^|,Oj , (12) 

and we have assumed a constant 

A In fc = In fcj+i — In ki . (13) 

The parameters pi then represent perturbations about 
the fiducial scale invariant model with (5<^ — 5.07 x 10^^. 
An arbitrary initial power spectrum is then represented 
as the linear interpolation in a log-log plot of points sam- 
pled at ki . If the initial power spectrum contains features 
below the sampling rate then the points represent aver- 
ages through the windows. This parameterization pre- 
serves the positive definiteness of the initial power spec- 
trum at the expense of making the initial power spectrum 
non-linear in the parameters. 

We can now define the infinitesimal response in the 
observable power spectra to the initial power parameters 
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This response function is a true transfer matrix for small 
perturbations pi from the fiducial model, 
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In Fig. (solid lines) we show this transfer matrix for 
three choices of Alnfc = 0.09, 0.05 and 0.02 and k, « 
km=Q ~ 0.13 Mpc~^. Note that the intrinsic resolution 
from the projection is only being reached at the finest 
binning for this high a k. 

The necessity of fine binning at high k rcfiects a trade- 
off in the choice of parameterization of the initial power 
spectrum. A logarithmic parameterization facilitates 
working with traditional approaches. Inflationary devia- 
tions from scale-invariance are usually quantified by the 
tilt n and its running a 

lnA?(fc) = ln<5? + (n-l)ln(-^) + ialn2(A)(i6) 
^ ^ fco 2 fco 
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FIG. 5: The power spectrum transfer matrix D^^ of a dis- 
crete perturbation in the initial power spectrum at in ki of 
width Alnfc. Here ki ~ km=e ~ 0.13 Mpc~^. The transfer is 
broadened by projection and lensing such that only features 
oi 5k/k > 0.05 can be resolved. 



Thus the traditional parameterization is a Taylor expan- 
sion in log- log space around the pivot point Infco- As 
a Taylor expansion, this parameterization can only be 
extended to a large dynamic range Aln(fc/fco) > 1 if 
a <^ {1 ~ n) <^ \ and can yield misleading results if ex- 
tended beyond its regime of applicability. Nevertheless, 
the effect of n and a can be easily constructed from the 
discrete parameters 



(n-l)ln(^) 



-a In ( — ), 
2 fco 



(17) 



to high accuracy with logarithmically spaced modes. 

The drawback of a logarithmic representation is that 
the geometric projection causes a smoothing that is 
roughly linear in i and hence the minimum A£/£ of fea- 
tures decreases with (.. However there is a second funda- 
mental limitation on features in the power spectrum that 
takes over at high ^: gravitational lensing 23]. 

Gravitational lensing deflects each line of sight accord- 
ing to the angular gradient in the projected potential 
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-$(x = nZ?,Z?,), (18) 



where the conformal lookback time Dc — D in the flat 
flducial model. $ is the Newtonian curvature and it is re- 
lated to the density perturbation in the comoving gauge 
by the Poisson equation. In the fiducial cosmology, an 
rms deflection of (t^ « 3' arises from the projected po- 
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tential 

al = Y.^{t+l)''-^Ct^ (19) 
i 

and gets its main contributions from £ ~ 50. In other 
words, the arcminute scale deflections have a coherence 
of several degrees. Also important are the fluctuations 
at the arcminute scale which cause a deflection of a few 
arcseconds. 

Following the linearized all-sky harmonic approach 
|24j . the change in the observed power spectrum is 

5cr' = ~w+i) -^^]^cr' (20) 

where se = and se — 2 and the mode coupling term 

.Mu,i, = ^[iiiii + 1) + e2{e2 + 1) - eii + 1)] (21) 
/ (24 + i)(2£2+"Ty f £ h i2\ 

""V in [s -s) ' 

a e + £1 + £2 = even; if odd. 

This mode coupling leads to a broadening of features 
in the angular power spectra in two ways. The main ef- 
fect is a large-angle warping of the temperature and po- 
larization maps that does not change the overall power 
in sub-degree scale fluctuations. Here the two terms in 
Eqn. H21|) mainly cancel and the linearized approximation 
retains validity beyond £^cr^/2 = 1 except for unphysical 
spectra with features of S£ ^ 1 ||23. However they in- 
duce a smoothing effect on features in the angular power 
spectrum and hence place a limit on resolving features in 
the initial curvature spectrum. In the absence of lensing, 
the resolution in harmonic space A£ is set by the funda- 
mental mode of the survey, i.e. for the full sky A£ — 1. 
In the presence of lensing, the incoherence of the deflec- 
tions in patches subtending more than a few degrees ef- 
fectively prevents their mosaicing to achieve higher reso- 
lution. This effect places a resolution limit of A£ 50 for 
the modes £ > 10'^ for which arcminute scale deflections 
are important. Thus a sampling rate of Alnfc — 0.05 is 
sufficient even on these scales. 

The second effect is from the smaller deflections with 
arcminute scale coherence. These induce additional ar- 
cminute scale fluctuations from the distortion of primary 
fluctuations on larger scales |25|. In the fiducial model, 
this secondary power becomes comparable to the primary 
power in temperature at ^soc ~ 3000 and essentially elim- 
inates the primary information on the initial power spec- 
trum beyond k > 4oc/-D* « 0.02 Mpc"^ Thus a sam- 
pling rate of Alnfc ~ 0.05 suffices for the whole observ- 
able regime. 

On the other hand the lensing potential power spec- 
trum Cf'^ is also a source of information about the curva- 
ture fluctuations at late times. Unfortunately due to the 
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FIG. 6: Individual errors on the discretized initial power spec- 
trum Pi with Alnfc = 0.05 and a cosmic variance limited 
measurement out to ^sec = 3000. With fixed cosmological 
parameters, temperature information alone (open triangles) 
recovers pi for modes near the acoustic peaks but not the 
troughs. Adding polarization information improves coverage 
to the whole acoustic range fc — 0.02 — 0.2 Mpc~^ (open 
squares). Adding cosmological parameters destroys individ- 
ual parameters constraints with just temperature information 
(filled triangles) which are largely recovered by adding polar- 
ization information (filled squares) . Much of the constraining 
power is hidden in the strongly anticorrelated covariance (not 
shown) . 

projection in Eqn. (|18|l . lensing has little ability to dis- 
tinguish sharp features in the initial power spectrum and 
even the interpretation of broadband power depends on 
the dark-energy and massive-neutrino dependent growth 
rate of structure '2^ . We ignore the information here by 
neglecting the sensitivity to the parameters coming from 
and the information in the B-mode polariza- 
tion. In this approximation, derivatives of the lensed 
CMB power spectra in Eqn. H14|l are simply given by the 
lensing of the power spectrum derivatives calculated from 
the transfer functions. In Fig. [S] (dashed lines) we show 
the effect of gravitational lensing on the power transfer. 
Now even features with A In fc = 0.05 are limited by lens- 
ing. We will hereafter adopt this spacing as the fiducial 
choice which leads to 200 parameters for the initial power 
spectrum. 

IV. STATISTICAL FORECASTS AND 
COSMOLOGICAL PARAMETERS 

A third limitation on the recovery of the initial power 
spectrum comes from degeneracy with cosmological pa- 
rameters such as the energy density of the particle con- 
stituents of the universe today. We employ the Fisher 
matrix technique to study the effect of degeneracies. 
Construction of the Fisher matrix will also be the first 
step in building statistics that are immune to cosmolog- 



7 



ical parameter uncertainty. 

We begin by appending the ric cosmological parameters 
to the initial power vector pi to form the joint parameter 
vector Pfj^, 11 = 1, . . . , Tife, . . . , Tifc + Tic- The Fisher matrix 
for a cosmic variance limited measurement out to ^ = .^sec 
is 

Ft.. = J2 ^^Tr[D,^C7iD,,C7i] , (22) 

where we have suppressed the XX' indices in a matrix 
notation for the derivatives and the power; here 

iBi,)xx- = ^f/' ^ 1^ (23) 

is the generalization of the transfer matrix Eqn. (|14|l . 
Given gravitational lensing constraints and other sources 
of foreground contamination, we will take ^scc = 3000. 

The inverse Fisher matrix approximates the covariance 
matrix of the parameters. In particular, the marginalized 
1-a errors on a given parameter cr{pfj_) are given by 

a'{p^) = C,,^{F-')^^. (24) 

The Fisher matrix also allows for a simple conversion 
between representations of the parameter set Pa and Pf_, 

where Pa can have a lower dimensionality, e.g. by going 
from the initial power parameterization of pi to n and a 
through Eqn. QTI) . 

In Fig. El we show the statistical errors on the ini- 
tial power spectrum parameters pi. Even in the absence 
of cosmological parameters {ric = 0) the constraints on 
the Pi are strongly modulated by the acoustic oscilla- 
tions with temperature information alone. The fc-modes 
associated with the acoustic troughs (half integral m) 
are poorly constrained while those associated with the 
peaks (integral m) are fairly well constrained. For low 
k, both the projection and cosmic variance severely lim- 
its the resolution in k leading to a strong anticorrelation 
between points. Hence only broadband power and not 
features can be recovered. For k > 0.2 Mpc~^ the cut off 
at €scc = 3000 limits the information. 

With the addition of the polarization and cross spectra 
information, the modulation of errors with the peaks re- 
verses. Now modes associated with acoustic troughs are 
better constrained than those associated with the peaks. 
The polarization information becomes even more impor- 
tant once the cosmological parameters are added back in. 
In Fig. El we take 6 additional parameters ^Ibh^, ^mh^, 
Ode, w, t and r, where r is the tensor-scalar ratio. With 
temperature information alone, the ability to measure an 
individual pi is essentially completely destroyed since the 
response to a cosmological parameter can a lway s be mim- 
icked by a particular superposition of pi's |27|. However 




0.1 0.2 0.3 



k (Mpc-') 

FIG. 7: First 5 principal components of the covariance matrix 
with 6 cosmological parameters, the amplitude In S( and tilt n 
marginalized. The modes are normalized to unit variance over 
dlnfc (discretized to Alnfe = 0.05). Well-constrained modes 
are localized in the acoustic regime k = 0.02-0.2 Mpc"^ and 
give zero weight to scale-free deviations. 

only six such combinations are lost and so with the addi- 
tion of polarization information with its complementary 
transfer function, the effect of marginalizing cosmological 
parameters is fairly small. In fact even with temperature 
information alone, there are modes that are very well 
constrained. The raw marginalized errors on pi are mis- 
leading in that they do not show the covariance. It is 
therefore useful to find a more efficient representation of 
the information that accounts for projection, lensing and 
cosmological parameters. 

V. PRINCIPAL COMPONENTS AND 
DEVIATED SPECTRUM 

Given the limitations imposed by projection, lensing 
and cosmology what kinds of deviations from a scale-free 
initial spectrum can the CMB best probe? To answer 
this question, we now consider a principal component 
representation of features. 

The information retained after including these three 
effects is quantified by the covariance matrix of the devi- 
ation parameters pi, the C'ij subblock of the joint covari- 
ance matrix C^^. Employing the Fisher matrix approxi- 
mation of the last section, we can determine the princi- 
pal component representation based on the orthonormal 
eigenvectors Sia of C'ij 

Q,=A-^Y.^^-''aSja, (26) 

a 
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a (eigennumber) 

FIG. 8: Eigenvalues and cumulative constraints on the run- 
ning of the tilt a represented in modes. For cosmic variance 
limited temperature and polarization measurements out to 
Isec = 3000, approximately 50 modes can be measured to 
percent level precision {oirna) < 0.01). Information on the 
running a saturates with more that ~ 10 modes. 

where A is a normalization parameter. For a fixed a, Sia 
specifies a linear combination of the deviation parameters 
Pi for a new representation of the spectrum 

Wg = A SiaPi , (27) 

i 

where the covariance matrix of the mode amplitudes is 
given by 

{mamb) = alSab ■ (28) 

This relationship can also be derived from Eqn. (|25|l . 
In other words, the eigenvectors form a new basis that 
is complete and yields uncorrelated orthogonal measure- 
ments with variance given by the eigenvalues. We choose 
the normalization A ~ (Aln/c)"^/^ so that a rria = 1 
corresponds to a deviation with unit variance when inte- 
grated over dlnfc, i.e. 

I d\nkSl^Y.^^a^^- (29) 

i 

Modes corresponding to combinations that are degener- 
ate with cosmological parameters or due to the projection 
and lensing have high eigenvalues and high variances and 
may be removed without loss of information. 

In Fig. [7| we show the first five eigenvectors of the 
covariance matrix including both temperature and polar- 
ization information. In this case we append In (5^ (or In A) 
and n to the list of cosmological parameters included in 
^IVI By doing so, the best constrained eigenvectors au- 
tomatically null out changes in the amplitude and tilt 
and so can be interpreted as the best modes to search 
for deviations from scale-free initial conditions. They 




k (Mpc-i) 

FIG. 9: Eigenvectors in the presence of instrumental noise 
(in ^K-arcmin) and beam (FWHM in arcmin). Parameters 
are chosen to roughly bracket the expectations of WMAP 
and Planck (dashed lines) for comparison with the noise-free 
case (solid lines). The eigenvectors are now confined to £- 
modes and hence fc-modes that are larger than the beam. In 
the highest noise, large beam case, polarization information 
is degraded and the avoidance of cosmological parameter de- 
generacies complicates the mode structure. 

also tend to null out any slowly varying component from 
residual secondary effects and foregrounds. As a conse- 
quence, their form is similar to Fourier modes confined to 
the well-measured acoustic regime 0.02 — 0.2 Mpc~^ with 
the constant and gradient term (in In k) removed. For the 
cosmic variance temperature and polarization spectra as- 
sumed in the Fisher matrix calculation, the mode errors 
are shown in Fig. |S1 (solid line) and remain below the 
percent level for ~ 50 modes. 

Note that truncation of the eigenmodes does not limit 
types of initial spectra that can be tested against the 
measurements nia since any spectrum may be projected 
onto the measurement basis. As an example, consider 
constraints on a constant running of the tilt a in the 
measurement basis. Each mode independently measures 
a with a variance 

alia) = Aj2sJ-ln\'^), (30) 

where the pivot fco is arbitrary since the Sia null out a 
constant deviation. Since each measurement is indepen- 
dent, the total variance is given by 

1 _ _J_ 

^'(") ~ h ^'(") ■ 

This cumulative sum is shown in Fig. |S| (dashed line). 
Information on a is efficiently captured in the first 10 
eigenmodes. 

The principal component technique can be made into 
a practical tool for searches for deviations from scale- free 



9 



initial curvature fluctuations. The best measured modes 
are fairly smooth and take forms that are independent of 
the discretization A In /c of the underlying representation. 
Hence they may be replaced by a continuous representa- 
tion for parameter estimation with an efficient integral 
code such as CMBFast or CAMB. The parameters 
can then be added in the usual way to a Monte Carlo 
Markov Chain. Deviations of the true model from the 
assumed fiducial model in constructing the eigenmodes 
and non-linearities in the parameter dependence will de- 
grade the statistical orthogonality of modes implied by 
Eqn. H28I) . These covariances will nonetheless be cap- 
tured in the likelihood exploration. The process may 
also be iterated if the assumed eigenvectors deviate sub- 
stantially from those of the best fit cosmology. 

Furthermore, although we have assumed a perfect ex- 
periment here to determine the limitations imposed by 
physical processes on feature recovery, the techniques are 
readily optimized for specific experiments. One simply 
replaces the power spectra in Eqn. (|22|l for the Fisher 
matrix with the total signal and noise power spectra of 
a given experiment. In Fig.EI we show a simple example 
where the experiment is parameterized by a white noise 
level in /uK-radian AT and the FWHM of a Gaussian 
beam 6. Here 



eei 



.AT, 



2^t{i+l)e^ ln2 



EE\ 



2(:^-)2gi'(«+l)e^/81n2 



T 



= 0. 



(32) 



We show the best eigenmode for two examples that 
roughly bracket the expectations for the completed 
WMAP and Planck missions. Note that the main ef- 
fect is to replace the cutoff imposed by tscc with that 
of the beam. However the smaller I range and the loss 
of polarization in the high noise example complicate the 
form of the mode because it adjusts to remain orthogonal 
to the cosmological parameters. 



VI. DISCUSSION 

Physical limitations imposed by the formation and evo- 
lution of CMB temperature and polarization anisotropy 
limit the ability to recover the initial power spectrum 



from even perfect measurements of angular power spec- 
tra. The initial power spectrum can be well measured 
in the acoustic regime k — 0.02 — 0.2 Mpc^^ but pro- 
jection and gravitational lensing prevents resolution of 
features smaller than Ak/k « 0.05. The complemen- 
tary acoustic polarization information greatly assists in 
filling in the coverage near the acoustic troughs and in 
immunizing the determination to cosmological parame- 
ter degeneracies. Furthermore polarization in principle 
offers a cleaner probe of long-wavelength power through 
reionization effects. Beyond purely two-point statistics 
in the temperature and polarization fields, the limita- 
tions of gravitational lensing can in principle be removed 
[28l l2^ . This can moderately improve the resolution of 
features but will only be possible with ultra-high preci- 
sion polarization measurements across nearly the whole 
sky. 

We have constructed the best statistics for measuring 
deviations from scale-free initial conditions given these 
physical limitations on power spectra. They are con- 
structed from the principal components of the covariance 
matrix. These modes resemble a Fourier expansion of the 
initial spectrum that is localized to the acoustic regime. 
Approximately 50 modes in this regime ultimately can 
be measured with percent level precision. These modes 
can provide many independent measurements of model- 
motivated parameters such as the running of the tilt as 
a consistency check that such simpler parameterizations 
are adequate. 

Although beyond the scope of this work, the princi- 
pal component technique can be developed into a practi- 
cal tool for the search for deviations from scale-free con- 
ditions with current experiments. These may also in- 
clude non-CMB sources with the appropriate addition of 
Fisher matrices. Here the physical limitations imposed 
on the modes would be supplemented by experimental 
limitations. The principal component technique provides 
a general, physically-motivated and complementary tool 
for the study of the initial power spectrum and its impli- 
cations for early universe physics. 
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